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Abstract 

Data from the Magsat spacecraft for November 1979 through 
April 1980 and from 91 magnetic observatories for 1978 through 
1982 are used to derive a spherical harmonic model of the earth's 
main magnetic field and its secular variation at epoch 1980.0. 
The model is called GS FC ( 1 2/83) * Constant coefficients are 
determined through degree and order 13 and secular variation 
coefficients through degree and order 10. The first degree 
external terms and cor respondi ng induced internal terms are given 
as a function of Dst. Preliminary modeling using separate data 
sets at dawn and dusk local time showed that the dusk data 

contains a substantial field contribution from the equatorial 

electrojet current. The final data set was therefore selected 
first from dawn data and then a. rented by dusk data to achieve a 
good geographic data distribution for each of three time periods: 
(1) November - December, 1979; (2) January - February, 1980; 

(3) March - April, 1980. A correction for the effects of the 
equatorial electrojet was applied to the dusk data utilized. The 
solution included calculation of fixed biases, or anomalies, for 
the observatory data. Although similar in many respects, 
GSFC Cl 2/835 differs from IGRF 1980 by 3.6 nT in the g^° term and 

shows a slightly negative B in the northern polar region as well 

as other differences in secular variation pattern. 
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Introduction 

NASA's Magsat spacecraft has provided the first truly global 
vector survey of the near-earth geomagnetic field. Data from 
early in the mission have previously been used to study the main 
Ci.e. originating in the earth's core) field in terms of an 
initial model (Langel et a l . , 1980), the spectral distribution of 
whe field (Langel & Estes, 1982) and as part of the data used in 
the 1960-1980 model of Langel et a l . (1982a) which was one of 
the candidates for the 1980 International Geomagnetic Reference 
Field (IGRF) . Since these early results the entire Magsat data 
set has been processed and a careful selection made of data 
suitable for modeling the main field. The present paper des- 
cribes our best estimate of that main field derived from the 
final Magsat data. 

Characteristics of the Magsat spacecraft, its instrumenta- 
tion and its operation have been published previously (Acuna et 
a l . 1973; Mobley et a l . , 1980; Lancaster et a l . , 1980; Langel et 

al., 1981, 1982b). The resulting vector data span the interval 

November 1979 through April 1980. After final attitude correc- 
tions and calibrations the data are estimated to be accurate to 
about 6nT (nanotesla) rss in each component and 2nT rss in field 
magnitude. This includes all known sources of error, including 
those of spacecraft position and attitude determination. The 
altitude range of the data is about 300-500 km. Orbital 
parameters were chosen so that Magsat was sun-synchronous at the 
terminator (i.e. the day-night boundary) so data are acquired 


uoen&m r 

OF POOR QU/Xii’j 


only at dusk and da* , Local times. 

The field is assumed to be curl-free and so represent ab l e by 
a potential function in the form of the usual spherical harmonic 
series: 

NMAX1 

v = r ° l 

n = 1 
NMAX2 n 

+ r y y (r/r_) n Cq,," 1 cos m cp + s m sin mcpH (cos 6-) 

o L. L o^n T n T n 

n = 1 m = o 


^(r n /r) n + Cg n in cos m cp + h n m sin mcp] P n m (cos$) 


o 

m * o 


(1) 

where r is the mean radius of the earth. taken to be 6371.2 km. 
o 

r, £ and cp are the standard spherical coordinates and p n m (cos -9- ) 
are the Schmidt quasi-normalised form of associated Legendre 
functions. The magnetic field is then given by 

B = - VV. (2) 

Theoretically (1) holds only if NMAXt and NMAX2 go to infinity 
and when the region of validity is source free. The measured 
internal 8 contains contributions both from the Earth's core and 
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from its crust and NMAX1 is chosen so that V represents fields 
from the core but not the crust, to our best estimation. Langel 
and Estes (1982) concluded that the core field dominates for n < 

13 and the crustal field for n >15 so, as in Langel et al. (1980 
and 1982a), we have chosen NMAX1 = 13. because Magsat passes 
through regions of "field aligned" currents, the source free 
assumption does not strictly hold. These currents have minimal 
effect on the field magnitude (Langel, 1974) so our procedure is 
to use the component data equatorward, but only the field 
magnitude poleward, of 50° geomagnetic latitude. 

The main contribution to the external portion of (1) comes 
from the equatorial ring current, with contributions also from 
magnetopause and magnetotail currents. Near the earth, fields 
from these sources tend to be aligned mainly along the dipole 
axis and so are well described when NMAX2 = 1. However, unlike the 
field from the core, the external fields vary considerably with 
both universal and local tiifie. The hourly Dst index is commonly 
taken to be an indicator of the relative change of these fields. 
Accordingly, Langel and Estes '( 1984 ) have investigated the 
relationship between q^° and Dst, which accounts for the 
universal time variation of It’s local time variation was 

investigated by carrying out the analysis separately for dawn and 
dusk giving the relationships: 


q^ 0 = 20.3 - 0.68 Dst 


Dusk: 


(nT) 


(3a) 
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Dawn: q., 0 « 18.62 - 0.63 Dst CnT) <3b) 

Similar relationships will hold at local timea inaccessible to 
Magsat. These are discussed using data from the POGO spacecraft 
in a later section of this paper. Because the external fields 
are time varying there are corresponding induced currents 

s 

within the earth which, in turn, contribute to the internal 
potential, e.g. g^ 0 . Accordingly g^° is expressed as a constant 
core field plus an induced internal field proportional to q.j 0 : 

Dusk: g^ 0 = -29987.7 + 0.24 q., 0 CnT) (4a) 

Dawn: g^ = -29992.3 + 0.29 q., 0 CnT) <4b) 

For greater validity in representation the present model includes 
relationships like (3) and (4), but omitting the local time 
dependance. 

Comgar^son of Dawn and Dusk results 

Division of the Magsat data into dawn and dusk subsets 
furnishes two independent data sets. Differences between 
models derived from these data sets should be an indication of 
the effects of local time asymmetries in external fields as well 
as in data quality and should also give some indication of the 
accuracy to which the spherical harmonic coefficients are deter- 
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mined. The data sets for dawn and dusk described by langel and 
Estes (1984 ) were modified slightly to improve geographic 
coverage and to eliminate a few data wi th -apparent external field 
effects. Their analysis was then repeated with nearly identical 
results. 

From equation (4) it is clear that there is a difference of 
about 5nT in the g^° term. Table 1 shows differences between 
coefficients of the amended models, up to degree/order five. 
Substantial differences in constant coefficients are present, 
especially for 9^°/ g 3°' 9 5°' ^3^ anc * ^ 5 ^ ' an< ^ t ^ ,e secular 
variation coefficients show even larger discrepancies. Differen- 
ces of this magnitude were totally unexpected and are too large 
to be due to data inaccuracy. A plot of the component differences 
indicates that the equatorial electrojet is the chief 
cause of the discrepancy. Plots of the dusk field minus the dawn 
field for the X (north) and Z (down) component are shown in 
Figures la and b. The dominant feature is a positive ridge in X 
following the geomagnetic equator with flanking positive Z to 
the north and negative Z to the south. This is precisely the 
result one would expect if the equatorial electrojet, which is 
eastward and below the spacecraft, were present at dusk and not 
at dawn. That the electrojet is present at dusk and not dawn is 
also found in our examination of individual and averaged 
passes in connection with studies of crustal fields and in the 
analysis of Roy (1983) who explicitly isolates and models the 
electro jet field. Maeda et al. (1982) and Roy (1983) also find 
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variations in the Y (east) component at dusk only, which they 
attribute to meridional currents. 

Figure 1 also shows a large dawn-dusk difference in the 
Antarctic (about 1?0° longitude, - 75° latitude) whose source, 
presumably ionospheric currents associated with the auroral 
belt, we are unable to identify with certainty. 

To model the field from the earth's core free from equa- 
torial electrojet effects requires either exclusive use of dawn 
data at the affected latitudes or correction of the dusk data. 
We chose the latter option because it made the final data 
distribution significantly better. The ideal correction would 
probably be a dynamic model of the equatorial electrojet. In the 
absence of any adequate modeL we resorted to the expedient of 
using the difference between the dawn and dusk potentials as the 
correction. Such a correction was applied to B and X between 
+ 20° geomagnetic latitude, to Y between + 15° geomagnetic 

latitude and to Z between + 50° geomagnetic latitude. It should 
be recognized that this correction removes a great deal of the 
i rcdependance of the dawn and dusk data sets. The model derived 
from this data set is referred to as the corrected dusk model. 

Columns 5 and 8 of Table 4 show the differences between the 
dawn and corrected dusk models and figures 2a and b show 
the new differences between calculated X and Z components. The 
coefficient differences no longer show isolated large values and 
are of amplitude in accord with the expected data accuracy. They 
should give some indication of the accuracy to which the 
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coefficients are determined. The difference plots no longer show 
evidence of the equatorial electrojet (as expected from the 
nature of the correction) but the higher latitude differences are 
still present. 

£i_nal Data Selection 

Because the model includes the secular variation it was 
necessary to obtain a good data distribution both in space 
and time. Accordingly/ after careful selection to ensure the data 
were from magnetically quiet periods (Langel and Estes, 

(1984) , the data were separated into three intervals: (1) 
November-December, 1979; (2) January-February, 1980; (3) 
March-April, 1980. For each period, and for dawn and dusk data 
separately, the data were collected into 5° x 5° equiangular bins 
over the globe. In regions with geomagnetic latitude equatorward 
of 50° where vector data with final attitude corrections were not 
available, field magnitude data were used from data whose 
attitude was known to a lesser accuracy. Residuals were computed 
from the GSFC(9/80) model (Langel et al, 1982) and out- 
lying points rejected. Within each bin data having 
outlying residuals relative to the mean for the bin were also 
re j ected. 

Data from each bin were then selected so as to obtain 
roughly the same number of points per equal area at all 
latitudes and to maintain an adequate spread of Dst values 
between + 22.5 nT to facilitate the external field analysis. 


Priority of selection went to the dawn data with dusk data, added 
to augment the spatial coverage. Dusk data within the 
"correction region" previously discussed were given a high 
weighting factor as were a few passes of lesser quality which 
were added to comp.lete the geographic coverage. Also a few 
vector points were selected from quiet polar cap data and 
added to the data set. Tablo 2 shows the distribution of selected 
data by component, local time and time period. 

Two versions of this data set were used: one with the 
correction applied to the dusk data and one with no correction. 

Preliminary modeling with the combined data set indicated 
that secular variation coefficients determined from 
Magsat data alone were significant only up to degree and order 
four. (This model is available upon request). Also, the coef- 
ficients so determined differ significantly enough from the 1980 
IGRF that we were unsure of their reliability. (See also Cain et 
al., 1983.) Accordingly we supplemented the Magsat data with the 
annual means from 91 observatories for the period 1978-1982. As 
described by Langel et al. (1982), the model included a solution 
for parameters representing the non-core field, or anomaly bias, 
at each observatory. Table 3 lists the observator i es used 
together with the anomaly bias from the solution and the 
difference between this anomaly bias and that from the GSFC C9/80) 
model. Data from each observatory were used after converting to 
X, Y and Z components in a geodetic coordinate system 
assuming an equatorial radius of 6398.165 km and reciprocal 


11 


flattening of 298.25 . The spatial distribution of the 
observatories Is displayed In Figure 4 . The observatory data 
were taken from the NOAA National Geophysical and Solar 
Terrestrial Data Center release 30 annual means data, with some 
updates for recently acquired data. This data set differs In 
some respects from the NOAA release 22 data set (updated by the 
authors) which was used for the GSFC(9/80) model In that the NOAA 
Data Center applied corrections to some observatories to account 
for discontinuities. Because of this sor.ie observatory biases in 
the present model differ from those for the GSFCC9/8C) model. 
The observator 1 es for which this is true are Alibag, Baker Lake, 
Bangui, Barrow, Bereznayki, Dikson and Kakioka. Only the 
observatories College, Fredericksburg, Guam, Honolulu, M'Bour, 
Maputo, Pamatai, and San Juan had-five annual means for the 
interval 1978.5 through 1982.5. Thirty-six observatories had 
four and the remainder three annual means. 

The Model 

The modeling method used is described in. Langel et al. 
(1982) . 

Both "corrected" and "uncorrected" models were derived but 
only the "corrected" version , called GSFCC12/83), will be 
presented here. Differences between the two were of the same 
magnitude as the difference between the dawn and corrected dusk 
models . 


The external field coefficients for GSFC(12/83) (not as a 
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function of local time since we have combined all data) are: 

q 1 0 = 18.4 - 0.63 Dst (nT) (5a) 
q 1 1 = -1.1 - 0.06 Dst (nT) (5b) 
s,j 1 » -3.3 *1* 0.17 Dst (nT) (5c) 

The corresponding g,j 0 internal coefficient/ including the effects 
of currents induced by the time-varying external field, is: 

g 1 0 = - 29991.6 + 0.270 q 1 0 (nT) (6) « 

Note that the Dst variation is determined only for the range - 20 
nT < Dst < 20 nT and may not hold outside this range. The 
corresponding induced contributions to g^ and h^ are regarded 
as ivjible and not computed. 

Coefficients for the internal field at 1980 are given in 
Table 4 along with the standard error determined in the fit and 
the difference between the dawn and corrected dusk model. We 
have contended before (Langel et al., 1982) that the standard 
errors of the coefficients as determined by the fitting process 
tend to under-estimate the coefficient inaccuracy. There are 
several possible reasons for this, such as deficiencies in the 
model validity, improper estimates of data correlation, and 
systematic effects of non-core fields. Prom Table 4, differences 
between dawn and dusk models are generally greater than the 


standard error of the coefficients. In the present case the two 
data sets should he** identical error characteristics/ except for 
the effects of unaccounted for fields from- ionospheric and 
magnetospher i c current systems* We expect/ then, that the coef- 
ficient standard errors reflect the internal consistency of the 
data whereas the differences between dusk and dawn models reflect 
the magnitude of the error from these external current systems. 
Table 5 summarizes the relative magnitude of the two sources. 
In general the estimated standard errors are too low by at least 
a factor of 2-3/ perhaps about a factor of 5 overall. The addi- 
tional error tends to occur at lower degree and/ from examination 
of Table 4/ within a particular degree at lower order. That is, 
the source is of long rather than short wavelength. Errors 
estimated on the basis of the coefficient differences may still 
be somewhat underestimated because of t ho toss of independence in 
the corrected dusk model. 


Secular variation terms up 

to degree and 

order 

ten 

are 

listed in Table 6 

together with 

the standard 

error 

of 

the 

coefficients. These 

ccef f i ci ents 

are determined 

i n a 

two 

step 


process wherin coefficients found insignificant in the first step 
are constrained to zero in the second step. (This procedure was 
not necessary for the main field model/ all coefficients were 
significant.) Differences between the degree and order four dawn 
and corrected dusk secular variation models were generally a 
factor of three or less larger than the corresponding standard 
errors. No such independant error estimate is available for the 
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final model incorporating the observatory data, but presumably 
the factor is similar. 

Table 7 gives the statistics of the various data sets 
relative to GSFC(12/83). The Magsat data used are before any 
selection required to achieve equal area distribution a/id so 
include more data than used in the model determination. As 
expected, the residuals are higher for the dusk data. The mean 
value of 4.3 nT in Y in the dusk data (and the 1.3 nT in the 
merged data) is due to the meridional current discovered by Maeda 
et al. (1982). The mean value of -4.6 nT in l in the dawn data 
(and the -3.7 nT in the merged data) is as yet unaccounted for. 
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Discussion 

Comparison of the G5FCC12/83) coefficients with those of 
earlier Magsat based models, including IGRF 1980, shows that the 
principal differences are the addition of the more accurate 
external field representation and of the secular variation 
determination. This is not to say that the small changes in 
constant coefficients are totally insignificant but rather that 
their effect is only seen in studies of small amplitude phenomena 
such as some of the details of crustal anomalies. An equally 
accurate independent data set is required to assess adequately 
the increased accuracy of this model. Because of the availa- 
bility of separate dawn and dusk data sets it was possible to 
give some idea of actual coefficient accuracy but this was 
impaired by the necessity of correcting the dusk data for the 
effects of the equatorial electrojet. Contour maps of the 
various components at 1980 are virtually indistinguishable from 
those of the IGRF (Peddie, 1982). 

For a model commensurate with the accuracy of the Magsat 
data inclusion of a model of external fields is a requirement. 
Table 8 gives a comparison of the g^° term from MGSTC6/80), 
GSFCC9/80) and IGRF 1980. Only MGSTC6/80) included a determina- 
tion of q^° and it did not include variation with Dst or take 
into account that some of g^ 0 could arise from induced currents.. 
The g^° difference between GSFC(9/8Q) (and IGRF 1980) and 
GSFCC9/80) is about -3.6 nT. If due totally to fields from 
induced currents this corresponds to a q^° of 13.3 nT and a Dst 



of 8.Q nT. The g^° difference between MGST (6/80) and 
GSFC(12/83) is much smaller. 

For some users of the model it may be desirable to take into 
account the local time variation of d-j°« Langel and Estes 
(1984) found that for data from the POGQ spacecraft: 

q^° = a+bDst (7) 

where a and b both varied with local time. The values of a at 
dawn and dusk were in disagreement with those from Magsat data 
but the values of b were in good agreement. It is suggested 
therefore that the expression for b from Langel and Estes be used 
in those applications where local time variations of q^° are 
important. That value is given by 

b = 0.680 + 0.134 sin (t + 13°), (8) 

where t is the local time in degrees r 

Recent examination of the Magsat data shows several small 
amplitude features which differ between local times but are 
relatively persistent for a particular local time. It is there- 
fore conceivable that a more detailed external field model is 
possible than that set forth here. 

To our knowledge the only previously published secular 
variation models for 1980 are GSFCC9/80) (Langel et al. 1982a), 
that of Barraclough et al. (1982), that of Peddie and Fabiano 
(1982), the IGRF 1980 (Peddie, 1982) and M061581 (Cain et al., 
1983). GSFC(9/8Q), while a good model for 1960-1980, is not 

suitable for prediction and will not be considered. IGRF 1980 is 
a combination of the models of °eddie and Fabiano (1982) and 
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Barraclough et at. (1982) and in our estimation is a good 

predictive model. Rather than considering each of Uhese three 
models we will discuss only IGRF 1980, 

Cain et a l . (1983) investigated the use of Mags at data, 

supplemented by Linear trends derived from quiet day midnight 
values at 19 observatories over the November 1979 to June 1980 

interval, in the determination of secular variation. They made a 

# 

comparison of B between M061581 and IGRF 1980. The 

largest discrepancy is in the northern polar region where M061581 
is strongly negative whereas IGRF 1980 is slightly positive. 
They attribute the difference to the effect of seasonal changes 
in external fields on the secular variation terms in M061581 and 
point out the possible importance of the effects of earth 
currents induced by changes in external current systems and of 
ionospheric currents below the satellite. 

B as derived from the degree/order four Magsat models (see 
earlier paragraph) also shows negative values at the northern 
pole. These are, however, greatly reduced in amplitude compared 
to M061581 (-20 to -40 nT/yr compared to -70 to -100 nT/yr), 

i ' 

presumably because the external field was more accurately 
accounted for in our model. Figure 3 shows plots of secular 
variation for the various components from GSFC(12/83). 
Comparison with plots from the IGRF (Peddie, 1982) shows a 

general agreement. There are, however, differences in some 

• 

details. The cell of large ' B in the Atlantic off of the 
Southeast Coast of the ll.S. has intensified by app rox i max e ly 10 
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nT/year and shifted slightly to the northeast. Also there has 
been a decrease in 3 over most of Europe of approximately 20 
nT/year. The pattern in the northern polar region closely 
resembles the structure of the degree/order four Magsat models 
discussed above. The zero nT/year contour extends to only the 
eastern coast of Greenland and turns back o«/er the Asian Arctic 
region, leaving the north polar region more negative (by 
approximately 20 nT/year) than indicated by IGRF 1980. Comparison 
of Z shows a decrease of app rox i ma t e ly 20 nT/year over central 
Asia, and a more complicated structure near the magnetic equator 
in the Pacific basin. The high positive cell off the coast of 
Equador has decreased in magnitude by 40 nT/year and moved to the 
east. To the ast of New Zealand, an intensification in Y of 
approximately 20 nT/year has appeared. The high negative X cell 
in the South Atlantic has intensified with an east-west 
signature, causing the region of low negative values off of the 
west coast of South America to disappear. 

The high accuracy of the constant coefficients determined 
from Magsat data does not extend to the secular variation model. 
Compared to IGRF 1980, GSFCC12/33) has the advantage of greater 
availability of observatory data spanning 1980 and of the availa- 
bility of a good (albeit short) temporal distribution of Magsat 
data. Statistics of the observatory data versus GSFCC12/83) are 
given in Table 7. Comparable statistics versus IGRF 1980 are 
found in Table 9. To make the two comparable, observatory 
"anomaly bias" values from GSFCC12/83) were used in both cases. 
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As expected the GSFC(12/83) model is a better representation of 
the data. 

Figures 5(a) - 5(c) show the yearly averages at a series of 
observatories together with the field predicted by GSFCC12/83) 
and IGRF 1980. These figures are fairly typical of the set of 
observatories used in the solution. On the whole, the secular 
variation trends predicted by the models are in close agreement. 
However it is clear that GSFC(12/83) benefits from the use of 

more recent data in defining the secular variation after 1980. 

# 

The Z variation of GSFCC12/83) at Muntinlupa (Figure 5(d)) shows 
the result of taking data over a short time interval when one of 
the annual means is not in agreement with the longer term trend. 
The secular variation solution is sensitive to the weights given 
the observatory data in the least squares adjustment. An initial 
model derived using the set of observatory weights from 
GSFC(9/8Q) produced first degree secular variation coefficients 

* o • 1 *1 

g^ = 23.3, g^ = 13.8, h^ = -20.9 and much less accurate 
secular variation trends at the observatories than the 
GSFC(12/83) model, which assigns a weight for each component at 
each observatory based on the standard deviation of a linear fit 
to the data over only the 1978-1982 interval. 

It is instructive to plot some of the spherical harmonic 
coefficients of secular variation from the various models avail- 
able in the literature as a function of time. In doing so we 
have tried to avoid selecting "predictive" models, i.e. models 
which are not based on data spanning the applicable interval. 


except for IGRF 1980. This eliminates several models for which 
the input data was graphically projected to a future epoch and it 
requires assigning a modified epoch to secular variation models 
based on data from a period of time earlier than the epoch of the 
associated main field model. Figure 6 contains plots of 9<|°/ 

- o • 1 *1 

^2 ' and ^2 ^ or t ^ ie “196 0-1980 time period from the models 

listed in Table 10. The dashed lines on the Figures are drawn by 
hand. The variation of g^° shows a change in slope (i.e. in 
secular acceleration) near 1968. The nature of the trend after 
1968 is unclear; there is some indication of a negative trend but 
if 6SFC(12/83) is correct the trend could either be constant or 
negative followed by positive. Since g 1 0 is the term most 
affected by solar cycle variations. perhaps there is a sunspot 
cycle variation (short dashes) superposed on the long term trend 
( long dashes) . 

A sharp change in trend is evident in g 2 at 1971/ after 

which the coefficient magnitude decreases sharply. A similar 

sharp change occurs in the h^ trend at 1971. These sharp 

changes in secular acceleration are further corroboration of the 

"1970 jerk" described, e.g., by Ducruix et al. (1980)/, Le Mouel 

et al. (1982), Malin et al. (1983) and Gubbins (1984). From 

Figure 6 we can obtain estimates of the secular acceleration and 

its change. For g 2 ° the pre- and post-1971 secular accelerations 

2 2 

are about -0.18 nT/yr and +0.80 nT/yr . Corresponding values 

12 2 

for h 2 are +1.0 nT/yr and -1.2 nT/yr . These values differ 

somewhat from those of Malin et al. (1983) but are in good agree- 
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merit with those of Gubbins (1984). The relation of the 1968 

0 A 

change in g<j to the "jerk'* is unclear. (Note / the values given 
by Malin et al. (1983) should be divided by four due to an 
algebraic error (Malin, private communi cation) which brings them 
into closer agreement.) 

Computation of "anomaly bias" values for observatory data 
was introduced by Langel et al. (1982) in an attempt to take into 
account the fields from crustal anomalies at observatory loca- 
tions together with any possible observer or instrumental bias. 
If no change in bias source occurs the bias should not change 
with time or from model to model. In practice, ability to 
determine accurately such a bias depends upon the adequacy of the 
model for the description of the true core field and its secular 
variation. For example suppose f(t) varies with time as: 

f(t) = b(t - t ) + c(t - t ) 2 (9) 

o o 

and we try to model f(t) with 

g(t) = a • + b « (t - t ) (10) 

o 

That is, f(t) has quadratic time dependance and no "bias" whereas 

the model is deficient in that it includes no quadratic time 

dependance but does include a bias. If f(t) is measured over the 

interval t - s/2 to t + s/2 and the square difference of f(t) 
o o 

and g(t) is minimized with respect to a* and b 1 over this 
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interval, then the solutions for a* and b 1 are b < « b and 
a» « (cs )/12. Thus, the inadequate model results in a false 
bias. 

Table 3 gives the difference in bias determination between 
G$FC<9/80) and GSFCC 12/83) for those observatories common to 
both. Note that the time periods over which the two determina- 
tions were made are disjoint as the GSFCC9/80) model -included 
observatory annual means only through 1977. The rms of these 
differences are given in Table 11 for each component together 
with the rms of the actual biases determined in the GSFC(9/80) 
model. Because the temporal variation over the four year time 
span of the GSFCC12/83) model is a great deal less complicated 
than that for the 20 year time span of the GSFCC9/8Q) model we 
expect that the temporal model of GSFC(12/83) should be more 
accurate. This leads us to conclude that the biases determined 
in the GSFCC12/83) model are probably more accurate. The dif- 
ferences are greater than we had expected and if they are an 
indication of model error rather than of a change at the observ- 
atory, then we are representing the local bias only to about 
15%. The large change in the Y bias at Bereznayki is related to a 
jump in the data which occurs at 1976.5, and the differences in 
bias values at Luanda Belas are related to the very rapid change 
in the data at 1981.5. These results may be caused by a site or 
instrument change or a change in procedures. The cause of these 
and other targe differences is under investigation. 
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Table 1; Difference between coefficients from the dawn spherical 
harmonic model and the dusk spherical harmonic model up 
to degree and order five 


n 

m 

An m 

A9 n 

u m 
Ah n 

, • fl! 
A9 n 

Ah n 

1 

0 

4.80 

0.0 

- 2.26 

0.0 

1 

1 

- 0.19 

-0.08 

2.65 

3.21 

2 

0 

0.87 

0,0 

2.50 

0.0 

2 

1 

0.20 

0.13 

-2.80 

0 . 77 

2 

2 

-0.06 

- 0.06 

1 .96 

- 0.68 

3 

a 

- 3.36 

0.0 

0.88 

Q.O 

3 

i 

-0.06 

- 1.08 

0.45 

3.27 

3 

2 

0.09 

- 0.03 

-1 .31 

0.60 

3 

3 

- 0.91 

0.25 

-1 .02 

- 0.10 

4 

a 

0.40 

0.0 

- 2.95 

0.0 

4 

i 

0.42 

0.36 

0.33 

- 2.62 

4 

2 

0.45 

Q . 1 3 

2.57 

-1 .88 

4 

3 

- 0.03 

- a. 21 

0.71 

0.15 

4 

4 

- 0.50 

0.22 

- 2.24 

1.77 

S 

0 

1.70 

0.0 

0.0 

0.0 

5 

1 

0. 25 

- 1.86 

o.a 

0.0 

3 

2 

- a . 45 

- 0.31 

0.0 

0.0 

5 

3 

0.17 

0.11 

0.0 

0.0 

5 

4 

0.06 

0.13 

0.0 

0.0 

5 

5 

- 0.42 

0.13 

o.a 

0.Q 


Table 2: 

Number of 

measu remen t s 

in final merged 

Data Sets 

Merged 
2 a _t a 


No v-Dec 

J an- F eb 

March-Ao 

X 

0 awn 

2320 

2591 

2596 


Dusk 

1833 

1609 

1620 

Y 

Dawn 

2251 

2457 

2419 


Dusk 

1902 

1775 

1 797 

Z 

Dawn 

2773 

3255 

3404 


Dusk 

1375 

948 

31 2 

a 

Dawn 

2975 

2832 

325 1 


Dusk 

2680 

2605 

2593 
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Table 3: Observatories Used in the Model 


Station 

Lat. Long. Alt. 

Anomaly bias (nT) 

Rias Oifferences 
with fiSFC(9/an) 



X y Z 

X Y Z 


ABISKO 

68.36 

18.82 

0.37 

19.3 

34.6 

29.7 



-28.3 

AOOIS ABABA 

9.03 

38.76 

2. 44 

307.3 

-27.6 

119.8 

17.8 

-42. 7 

ALERT 

32.30 

-62.30 

0.06 

-6. 1 

31.2 

-197.3 

-11.3 

-3. 4 

-34.3 

AC I BAG 

13.64 

72.37 

0.0 

-207. 1 

437. 9 

394. 4 

-126. 9 

-5.8 

—7 C . 6 

ALMA ATA 

43.23 

76.92 

1.30 

148.6 

37.8 

-183,8 

MO -J 
0m Jtm 0 f 

-4. 4 

-31.7 

ALMERIA 

36.33 

-2.46 

0.06 

-14.3 

17,3 

3.4 

-8.3 

8.9 

-16.3 

AMATSIA 

31.33 

34.92 

0.0 

103.0 

38.2 

273. 1 



US. 1 

ANNAMALA I NAGAE 

11.37 

79.68 

0.0 

178.8 

-90.3 

-39.3 

-40.7 

-9.9 

AQUILA 

42.38 

13.32 

0. 62 

-1.3 

37.3 

-3, 1 

7. 1 

-4. 7 

-12.6 

ARGENTINE ISCNO 

-63.24 

-64.26 

0.0 

68.3 

-74.6 

486.7 

-1.0 

8.9 

-26. 1 

BAKER CAKE 

64.33 

-96.03 

0.04 

161.3 

-47. 1 

-92. 4 

-37.9 

62.6 

21.9 

BANGUI 

4.44 

18.37 

0.38 

-148.3 

-26. 1 

222.9 

-29.0 

-62.7 

113.6 

BARROW 

71,32- 

136.62 

0.0 

21.6 

-69. 7 

-33. 4 

-2.3 

-3. 6 

-23. 4 

BELSK 

31.34 

20.79 

0. 18 

107.7 

132.3 

298.7 

11.2 

-21.2 

-13.3 

bereznayki 

49. 82 

73.08 

0. 0 

-379.6 

-193.2 

310.6 

44.2 

-207. <= 

-28. 3 

BJORNOYA 

74.30 

19.20 

0.07 

-110.0 

37.4 

18.0 

12. 1 

10.3 

-10. C 

BOULDER 

40. 14-103.24 

1.63 

-14.2 

30.6 

-166.8 

38.3 

-2.9 

-6.3 

CAMBRIDGE BAY 

69. 20- 

103.00 

0.02 

102.7 

-96.9 

116. 1 




CASEY 

-66.23 

1 10.33 

0.0 

906.2 

-307.9 

-669.4 




CHANBON FORET 

4G.02 

2.26 

0. 1* 

-71.8 

-21.7 

103.3 

7,9 

2- 9 

-13.2 

COrfIBRA 

40.22 

-8.42 

0.09 

21.3 

-3.3 

3.3 

2o. 7 

17.3 

-42.7 

COLLEGE 

64. 36- 

147. 84 

0.09 

-18.0 

-33.3 

-106.8 

-0. 3 

3.2 

0. 3 

DIKSON 

73.34 

80.36 

0.01 

-83.0 

-133.7 

-262.8 

19.5 

13.9 

-14.2 

DOURBES 

30. 10 

4.60 

0. 20 

4. 0 

-16.6 

81.8 

11.7 

0, 2 

-16. 6 

DUnONT DURVILLE 

-66. 66 

140.01 

0.03 

-130. 1 

-403.8-2861.2 

4 J . 9 

-13.3 

-8. 3 

OUSHETI 

42.09 

44.71 

0. 98 

-221.3 

6.3 

-114.3 

16. 9 

1.7 

-31.3 

ESKDALEMUIR 

33.32 

-3.20 

0.24 

4. 7 

-49.4 

-32.6 




FORT CHURCHILL 

38. 77 

-94. 10 

0.04 

-116.3 

31 

-269.0 

23. 7 

-13.2 

-13.2 

FREDERICKSBURG 

38.21 

-77.37 

0.06 

36.9 

-63.9 

133.4 

23.8 

5.3 

14. « 

FUR3TNFELDBRUCK 

48. 16 

11.28 

0.37 

-24.8 

-8.3 

12. 1 

6.4 

-2. 4 

-2.3 

GN ANGARA 

-31.78 

1 13.93 

0.06 

-3.4 

-126.7 

142.3 

47,6 

t. 1 

19. S 

OORNOTAYEZHNAYA 

43.68 

132. 17 

0.30 

-8.9 

-11.2 

-62.2 

4. 4 

12.3 

-24. 4 

GREAT WHALE R 

33.27 

-77.78 

0. 02 

261.0 

102.0 

-79.7 

2. 7 

-11.5 

-20. 8 

GUAM 

13.38 

144.87 

0. 13 

133.6 

86.2 

76.6 

-41.3 

-18. 1 

79.3 

HARTLAND 

30. 99 

-4. 48 

0. 09 

-42. 6 

11.6 

62.2 

12.0 

4.6 

-2. 0 

HER MANUS 

-34. 43 

19.23 

0.02 

18.2 

11.3 

24.3 

33.8 

-12.7 

30.2 

HGNOLULU 

21.32- 

138.00 

0. 0 

-167.3 

80. t 

-340.9 

6. 9 

-32.8 

-36. 3 

HUANCAYO 

-12.04 

-73.34 

3.31 

69. 4 

36. 1 

17.3 

-20.9 

-12.4 

-229.3 

HYDERABAD 

17. 41 

78.33 

0. 49 

303.3 

20.8 

476.3 

-46.9 

-88.3 

19.3 

KAKIOKA 

36.23 

140. 19 

0.02 

-11.8 

10.4 

-86.6 

-3.0 

7.3 

— 2 n , 3 

KANOYA 

31.42 

130.88 

0. 10 

-13. 1 

33.3 

-37.3 




KERGUELEN 

-49. 33 

70. 20 

0.04 

223.0 

189.8 

643. 4 

43. 3 

-16.3 

10. 4 

KODA IKANAL 

10. 23 

77.46 

2.32 

-348.8 

287.9 

-32. 1 




LEIRVOGUR 

64, )8 

-21.70 

0.0 

-283. 1 

601.3 

-488. 3 

14. 3 

6. 4 

-3.o 

LERWICK 

60. 13 

-t. 18 

0.08 

-127.2 

169.9 

33.3 

10. 1 

0.6 

-2. 2 

LOVO 

39.33 

17.83 

0.02 

44 . 8 

-4.8 

-3.3 

°. 6 

-3.7 

-11. 0 

LUANDA BELAS 

-8.92 

13. 17 

0.03 

290.3 

-30. 1 

204.8 

A7.3 

-23. 7 

134.3 

LUMPING 

23.00 

121. 17 

0. 10 

6.6 

43. 7 

40.0 

-9. 4 

30.3 

-39. * 

LVOV 

49.90 

23. 73 

0.20 

140.2 

121.4 

146.3 

11.3 

-3.9 

-7.4 

fl BOUR 

14.39 

-16. 96 

0.0 

119. 9 

46.8 

33.3 

-23. 1 

95.2 

-23. 3 

MACQUARIE tSLNO 

-34.30 

138. 93 

0.0 

274.3 

2.4 

292.3 

38. 0 

9.7 

-16. 7 

MAPUTO 

-23. 92 

32.38 

0.03 

349.7 

34.0 

-137.3 




MAUSON 

-67. 60 

62.38 

0.0 

23.3 

21. 9 

132. 0 

6. 6 

-3. 1 

-13. 7 

MEANOOK 

34 . 62- 

113. 33 

0.6/ 

1 18.4 

12.6 

-142.0 

30.5 

4,3 

18.3 

memambetsu 

43. 91 

144. l? 

0. 03 

-243. 4 

141.9 

63. 1 

-4. 1 

S. 8 

-12. 1 

MIRNYY 

-66. 33 

93. 02 

0.02 

-102.9 

48.3 

-430. 1 

32.2 

-l.« 

-2i.o 

MOULD BAY 

76.20- 

119.40 

0. 13 

-20.9 

-1.6 

-63.3 

1.3 

—6. 6 

-24.6 

MUNTINLUPA 

14.38 

121 .01 

0.06 

- 44 . J 

-44.3 

33. 4 

20. 1 

139.4 

16. 3 

NAGYCENK 

47 . 63 

16. 72 

0. 15 

13.0 

-3. 6 

-31.0 




NEWPORT 

48.26- 

117. 12 

0.77 

-38.8 

111, 1 

-123.0 

32. 5 

0.9 

6.2 

NIEMEGK 

32.07 

12.68 

0.07 

-31.9 

-1.4 

-82. 4 

9.4 

-1.5 

-14.4 

NURMI JAR VI 

60.31 

24. 63 

0. 11 

281.8 

-106. 1 

93. 3 

10. 7 

-7.3 

-l 1 . 

OTTAWA 

43. 40 

-73.53 

0.76 

133.9 

-146.3 

137.3 

29 . 2 

-3. 3 

-27. 4 
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Table 3: Observatories Used In the Model 


-v I 


Station 

Lat. Long. 

Alt. 

Anoma 

X 

ly bias 
y 

(nT) 

Z 

Bias Differences 
with GSFC(9/R0) 

X Y Z 

PA«A7AI 

-17.37-149.37 

0.09 

-664, 3 

-732.2 

-127. 7 

11.3 

-40.7 

32.0 

PLESHENITZI 

34.30 27.88 

0.20 

273.8 

171.7 

-139.8 

6. 4 

23.4 

-20. 2 

PORT MORESBY 

-9.41 147.13 

0.08 

11.3 

49.3 

273.3 

31.7 

-10.3 

109.9 

PORT-ALFRED 

-46.43 31.37 

0.0 

-696.4 

1 137. 7 

133.9 




QUETTA 

30.19 66.93 

1.73 

-8. 4 

83.2 

-34. 7 

2_. o 

32.3 

-0. 1 

RESOLUTE BAY 

74.70 -94.90 

0.03 

43.2 

29.0 

36. 1 

3.3 

-1. 1 

-13.9 

SABHAUALA 

30.36 77,30 

0.49 

-19.6 

-63. 9 

13.0 

-4. 4 

-36. 1 

-48.0 

SAN JUAN 

18. 11 -66. 13 

0.40 

-30.9 

l T y . 4 

133.8 

88. 4 

-4.6 

-29.3 

SANAE 

-70.30 -2.37 

0.03 

-24. 7 

3. 7 

14.3 

26.7 

27.3 

-14, 1 

SITKA 

37.06-133.32 

0.02 

-3.4 

-19.8 

-67.2 

6.2 

-3.3 

13. 7 

SODANKYLA 

67.37 26.63 

0.17 

-139.4 

-106.4 

-339.9 

3.2 

-2. 1 

-11.0 

ST JOHN S 

47.60 -32.68 

0.0 

37.2 

22.6 

-6.0 

48.2 

Q 

-3. 9 

SURLAR I 

44.68 26.23 

0.08 

10.3 

-32. 1 

-69. 1 




THULE II 

77.48 -69. 17 

0.03 

-63.4 

103.3 

3.0 

-7.3 

-1.4 

-S.3 

TIHANY 

46.90 17.39 

0. 19 

-26.3 

6.3 

-43.3 

-10.2 

•>-Tv A 

*. • ** 

-16 l 

TOLEDO 

39.38 -4,03 

0.30 

9.7 

4. 4 

-3.0 




TRIVANDRUM 

8.48 76.93 

0.29 

276.7 

204. 1 

210.3 

21.7 

8.9 

137.3 

TROMSO 

69.66 18.93 

0. 11 

l 10. 1 

-407.3 

l 17. I 

34.6 

-3.3 

-30.3 

TSUMfB 

-19,22 17.70 

0.08 

33.0 

-68.0 

l 13.6 

38.9 

-8. 0 

32.6 

TUCSON 

32.23-110.83 

0.76 

-*1 . 1 

-39.0 

133.2 

33. 0 

3. 1 

-33.9 

VALENT I A 

31.93 -10.23 

0.0 

134.2 

-46.0 

19. 3 

13. 9 

10. 1 

-6. 1 

VANNOVSKAYA 

77.93 38.11 

0.37 

170. 1 


64,6 

-l . 4 

3. 0 

3.8 

VASSOURAS 

-22.40 -43.63 

0.43 

78.0 

-60. 3 

-33.0 

33.3 

66. 0 

36.9 

VICTORIA 

48.32-123.42 

0. 19 

31.8 

1.7 

-327.3 

49. 4 


0. 7 

UHITESHELL 

49.73 -93.23 

0.0 

132.3 

-242.3 

-230. 1 




WINGST 

33.74 9.07 

0.04 

33.3 

43. 1 

-66. 3 

.10.9 

-3. 0 

-l 1 . 3 

WITTEVEEN 

32.31 6.67 

0.01 

23.3 

-0.2 

-78. 2 




YELLOW-KNIFE 

62. 40-114.30 

0. 18 

401.2 

-216.0 

129. 7 
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aKoefficients for Epoch 198Q, Units are ?T. mean radius of the 
earth is 6371 „2 km 


n 

(til 

«„■ 

Standard 

error 

Dawn-Dusk 
Di f f erence 

h m 

h n 

standard 

error 

Dawn-Dusk 
D i f f erence 

l 

0 

•- 2999 1 . 6 

0 . 05 

-0 . 30 




1 

1 

-1956.0 

0.03 

-0 . 4 1 

5603 . 9 

0 . 08 

-0 . SO 

’T> 

0 

-1996.7 

0 . 05 

-0 . 33 




O 

1 

3027.3 

0.05 

-0 . 3 1 

-2129.3 

0.05 

-0 . 1 0 


»•» 

1662.7 

0.07 

-0 . 03 

-199.7 

0 . 07 

0.08 

3 

0 

1231.4 

0.05 

0.33 




■*i 

•Z* 

1 

-2 ISO. 5 

0 . 06 

0. 16 

-335.5 

0 . 06 

0 . 44 


2 

1250.9 

0.05 

0.14 

270.9 

0 . 05 

-0 . 06 

o 

3 

332 . 9 

0 . 06 

0. 15 

-252 • 2 

0. 06 

-0 . 09 

4 

0 

937.6 

0.05 

-0.20 




4 

1 

732.3 

0.05 

-0. 17 

212.2 

0 . 05 

-0 . 37 

4 

2 

397 . 9 

0 . 05 

-0 . 07 

-256 . 7 

0 . 05 

-0.07 

4 

3 

-413.9 

0.05 

0 . 1 1 

' 52.9 

0 . 05 

0 . 1 2 

4 

4 

193.9 

0.06 

0. 13 

-297.2 

0 . 06 

-0. 07 

ej 

0 

-217.9 

0.05 

0.03 




5 

1 

357. 2 

0.05 

0 . 09 

46.0 

0 . 05 

0 . 67 

5 

2 

261.0 

0.05 

0.40 

149.3 

0 . 05 

0 . 1 4 

5 

3 

-74.3 

0 . 05 

0.05 

-150.6 

0 . 05 

-0 . 09 

5 

4 

-161.9 

0 . 05 

0.01 

-77.7 

0 . 05 

-0.04 

S 

5 

-43.0 

0.06 

0.16 

92. 1 

0 . 06 

-0. 06 

6 

0 

43.0 

0. 05 

0 . 23 




6 

1 

65.5 

0 . 05 

0 . 29 

— 14.3 

0.05 

-0.34 

6 


41.9 

0.05 

-0 . 04 

93. 1 

0 . 05 

-0 . 26 

6 

3 

-192. 1 

0.05 

-0 . 06 

70.6 

0 . 05 

0. 10 

& 

4 

3.6 

0 . 05 

0.02 

-43. 1 

0.05 

-0 . 03 

6 

5 

13.3 

0.05 

-0.01 

—2 . 2 

0 . 05 

-0 . 06 

6 

6 

-107.7 

0.06 

-0 . 02 

17.2 

0 . 06 

-0 . 03 

7 

0 

71.9 

0.05 

-0 . 53 




/ 

1 

-59. 2 

0.05 

-0. 14 

-32.5 

0 . 05 

0 ■ / 

7 

2 

1.7 

0 . 05 

0 . 54 

-27 . 4 

0 . 05 

0.53 

7 

3 

20.7 

0.05 

0 . 06 

-4.9 

0 . 05 

-0 . 30 

7 

4 

-12.4 

0 . 05 

-0. 16 

16.3 

0 . 05 

-0 . 05 

7 

5 

0.6 

0.05 

0 . 03 

17.3 

0.05 

-0 . 05 

7 

6 

10.6 

0.05 

0 . 0 1 

-23 . 0 

0.05 

-0. 03 

— * 

/ 

7 

-1.7 

0.06 

0 . 03 

-9.7 

0 . 05 

-0 . 02 

P. 

*0 

13.5 

0.05 

-0 . 06 




S 

1 

6.5 

0.05 

0.26 

6 . 3 

0. 05 

—O ■ 1 4 

8 


-0.4 

0.05 

-0 . 0 1 

-17.7 

0 . 05 

-0 . 3 ■> 

3 

3 

-11.0 

0 . 05 

-0. 03 

4.2 

0 . 05 

-0 . O'. * 

3 

4 

-6.9 

0.05 

-0 . 04 

— O 

0 . 05 

0 . 00 

8 

5 

4.2 

0.05 

-0.04 

9. 1 

0.05 

-0 . 02 

3 

6 

2.7 

0.05 

-0.03 

16. 1 

0. 05 

-0 . 00 

3 

7 

6.0 

0.05 

0.04 

-13.2 

0.05 

-0 . 0 1 

s 

3 

-1.4 

0 . 06 

-0 . (3 1 

-14.6 

0.06 

-0 . 05 

9 

0 

5.3 

0.05 

0 . 23 




9 

1 

10.4 

0 . 05 

0.18 

-20 . 8 

0 . 05 

0 . 06 

9 

P 

1.4 

0 . 05 

0. 27 

15.5 

0 . 05 • 

0 . 23 

9 

3 

-12.3 

0.05 

0. 13 

3.7 

0.05 

-0.25 

9 

4 

9.4 

0.05 

-0. 10 

—5. 3 

0.05 

-0. 14 

9 

5 

-3.4 

0.05 

-0 . 00 

-6.3 

0.05 

0.09 

9 

6 

-1.2 

0.05 

0.01 

9.0 

0.05 

0.07 

9 

7 

6.7 

0.05 

0.09 

9.6 

0.05 

0.07 

9 

3 

1.5 

0.05 

0.03 

-6.0 

0.05 

-0.03 

? 

o 

-5.0 

0.06 

-0. 02 

2.0 

0.06 

-0.03 


31 


Table 4 Cont'd 


ORIGINAL fW2 rJ 
OF POOR QUALITY 


n 

« 

- * 

Standard 

Dawn-Dusk 

h n" 

standard 

Oawn-Ousk 



*n 

error 

Difference 


error 

D 1 f f erence 

10 

0 

-3. 5 

0.05 

-0.23 




to 

1 

-4.0 

0.05 

0,07 

1 .2 

0. 05 

0.24 

10 


•7* O 

Jim ■ 

0.05 

0.06 

0 . 5 

0. 05 

-0. 17 

10 

3 

-5. 4 

0 . 05 

-<5.14 

2.7 

0.05 

-0 . 04 

10 

4 

-2.0 

0. 05 

0 . 02 

5 . 8 

0 . 05 

-0 . 0 1 

10 

3 

4.6 

0. 05 

-0.01 

-4.2 

0 . 05 

-0 . 05 

10 

6 

3. 1 

0.05 

-0.03 

-0 . 4 

0. 05 

0 . 02 

10 

7 

0.9 

0.05 

0.02 

-1.3 

0 . 05 

0 . 08 

10 

8 

2,0 

0.05 

-0.08 

3.6 

0. 05 

-0. 04 

10 

9 

2 . 3 

0 . 05 

0.02 

-0 • 5 

0.05 

0.04 

10 

10 

-0.3 

0 . 00 

-0.01 

-6. 1 

0 . 06 

-0 . 03 

1 1 

0 

2.5 

0 . 05 

0. 14 




1 1 

1 

-1 . 1 

0.05 

0. 10 

0 . 6 

0.05 

-0 . 1 2 

1 1 

*“» 

«!• 

-1 . 7 

0.05 

0 . 05 

1.7 

0 . 05 

0.15 

11 

•_* 

£.2 

0.05 

0. 15 

-1.3 

0. 05 

-0 . 02 

11 

. 4 

0. 1 

0 . 00 

0 . 04 

-3. 1 

0 . 05 

-0 . 1 9 

11 

5 

-0.6 

0 . 05 

0 . 05 

0 . 8 

0. 05 

0 . 03 

1 1 

6 

-0 . 3 

0 . 05 

0.01 

-0 . <5 

0 . 05 

0 . 02 

1 1 

7 

1 . 6 

0 . 05 

0. 10 

—2. 5 

0 . 05 

0.02 

1 1 

3 

1.7 

0 . 05 

0 . 08 

-0.4 

0.05 

0 , 0 1 

1 1 

9 

-0 . 3 

0. 05 

-0.00 

-1.7 

0 . 05 

-0 . 09 

1 1 

10 

2.0 

0.05 

-0.01 

-1.5 

0 . 05 

0 . 0 1 

1 1 

11 

3.4 

0.06 

-0.02 

0 . 8 

0.06 

-0 . 0 1 

12 

0 

-1.7 

0 . 05 

-0. 12 




12 

1 

-0 . 1 

0 . 05 

-0. 12 x 

0.2 

0.05 

0 . 20 

12 

'7* 

-0 . 3 

0 . 05 

0.22 

0.8 

0 . 05 

-0 . 05 

12 

•J 

-0.2 

0.05 

-0 . 05 

2.5 

0. 05 

-0 . 02 

12 

4 

0.7 

0.05 

-0 . 02 

-1.5 

0 . 05 

-0 . 05 

12 

5 

0 . 6 

0.05 

-0.02 

0 . 4 

0 . 05 

-0 . 03 

12 

6 

-0. 4 

0. 05 

0 . 00 

0.2 

0 . 05 

0 . 05 

12 

7 

-0. 2 

0. 05 

-0 . 00 

-0 . 2 

0 . 05 

0 . 08 

12 

s 

0.3 

0 . 05 

-0.03 

0 . 0 

0 . 00 

0 . 02 

12 

9 

-0.5 

0.05 

<5 . 03 

-0.0 

0 • 00 

0 . 02 

12 

10 

0. 1 

0.00 

-0.01 

-1.2 

0 . 05 

— 0 . 03 

12 

1 1 

0.6 

0.05 

0.02 

0.4 

0 . 05 

-0.04 

12 

12 

-0. 1 

0.00 

0.02 

0. 5 

0 . 06 

-0 . 06 

13 

0 

-0 . 2 

0.00 

0. 10 




13 

1 

-0 . 6 

0 . 05 

<5 . 03 

-0.5 

0 . 05 

-0. 03 

13 

•? 

0. 4 

0. 05 

0 . 08 

0 . 3 

0 . 05 

0. 10 

13 

3 

-0 . 9 

0 . 05 

0. 14 

1.5 

0 . 05 

-0 . 0 1 

13 

4 

-0.2 

0 . 00 

0 . 06 

-0 . 2 

0.05 

-0 . 1 1 

13 

5 

1 . 1 

0. 05 

0 . 05 

-0 . 5 

0 . 05 

-0 . 03 

13 

6 

-0.4 

0. 05 

-0.05 

-0 . 2 

0 . 00 

0 . 08 

13 

7 

0.4 

0 . 05 

0.09 

0 . 8 

0 . 05 

0 . 06 

13 

8 

-0.5 

0 . 05 

0.06 

0. 1 

0 . 00 

0 . 04 

13 

9 

1 0.2 

0 . 05 

0 . 02 

0 . 9 

0 . 05 

-0 . <1)8 

13 

10 

-0. 1 

0 . 05 

-0. 03 

0. 1 

0. 00 

0 . <1)3 

13 

11 

0. 3 

0 . 05 

0 . 0 1 

-0 . 1 

0 . 00 

0 . 09 

13 

12 

-0. 1 

0 . 00 

0 . 03 

0 . 3 

0. 00 

0. 04 

13 

13 

0. 4 

0 . 06 

0 . 03 

— 0 . 3 

0 . 06 

— 0 . <1)4 
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Rms of 

Rms of 

Degree 

standard 

Dawn-Dusk 


error CnT) 

Difference CnT) 

1 

0.07 

0. 69 

2 

0. 06 

0. 40 

3 

0 . 06 

0 .' 23 

4 

0 . 05 • 

0. 17 

5 

0.05 

0 . 25 

6 

0 . 05 

0. 16 

7 

0 . 05 

0 . 34 

a 

0 . 05 

0. 1 1 

9 

0.05 

0 . 1 5 

10 

0 . 05 

0, iO 

11 

0 , 05 

0.03 

12 

0 . 05 

0. 03 

13 

0.04 

0 . 06 


Ratio 


9 . 9 
6.7 

•-« • C» 

3.4 

5.0 
3.2 
6 • 3 

3 . 0 

2 . 0 
1.6 
1.6 

1 . 5 



t 


i 


Table 6 : The GSFCC12/33) secular variation model 
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a) Coef f i ci ents for Epoch 1930/ Units are nT/yr mean radius of the 
earth is 6371 . 2 km 


n 

m 

5 * 

Standard 

Dawn-Dusk 

u » m 
h. 

s t anda rd 

Dawn-Ous k 



n 

error 

Di f f erence 

n 

error 

Di f f erence 

1 

0 

26.5 

0. 17 

1 . 66 




1 

1 

11.9 

0 . 23 

-0. S3 

-16. 0 

0 . 20 

-1.37 

2 

0 

-17.S 

0. 15 

- 1 . 03 




? 

1 

3.6 

0. 16 

1 . 00 

-15.0 

0. 19 

0 . 29 


*"> 

3. 1 

0 . 22 

-1. 15 

-24. 2 

0 . 22 

0. 10 

3 

0 

1.2 

0 . 1 3 

-0 . 03 





1 

— * Z * 9 6 

0. IS 

-0.49 

o C; 

0. 15 

-2.44 



-1. 1 

0. 1 7 

0. 10 

2.5 

0. 15 

0 . 02 


3 

1 . s 

0 . 20 

0 . 6 1 

-4. 1 

0 . 22 

0 . 03 

4 

0 

-0. 7 

0. 14 

0. 92 




4 

1 

-2. 1 

0. 13 

0 . 29 

2 m 9 

0. 14 

1 . 49 

4 

'“i 

6T C» 

• o 

0. 15 

- 1 . 93 

2. 0 

0. 14 

0. 77 

4 

i a i 

-2. 4 

0. IS 

-0 . 33 

2.0 

0. 16 

0 . 29 

4 

4 

__ o 

*1 ■ o 

0.25 

0. 69 

-2. 4 

0. 19 

0 . 07 

5 

0 

0.0 

0 . 00 

0 . 00 




5 

1 

-0.0 

0 . 00 

0 . 00 

1.4 

0. 15 

0 . 00 

3 

O 

-1.3 

0. 12 

0 . 00 

0 . 5 

0. 13 

0 . 00 

5 

•r» 

C* 

-4. 1 

0. 13 

0 . 00 

-1.0 

0. 14 

0 . 00 

5 

4 

-1.3 

0. 17 

0.00 

0 . 6 

0. 15 

(!) . 00 

5 

5 

0 . 0 

0 . 00 

0 . 00 

0. 1 

0. 17 

0 . (DO 

6 

0 

1.6 

0. 12 

0 . 00 




6 

1 

0. 4 

0 . 09 

0 . 00 

0.4 

0. 10 

0 . 00 

& 

O 

1.0 

0.11 

0 . 00 

-0.7 

0. 10 

0 . 00 

6 

3 

1.3 

0 . 09 

0.00 

0 . 0 

0 . 00 

0 00 

6 

4 

0 . 6 

0. 15 

0. 00 

-0. 4 

0. 12 

(!) . 00 

6 

5 

0.3 

0. 10 

0 . 00 

0.4 

0. 13 

0. 00 

6 


1.6 

0. 14 

0 . 00 

1.6 

0. 16 

0 . 00 

— «r 

/ 

0 

1 . 4 

0. 10 

0.00 




7 

1 

-0. 7 

0. 03 

0. 00 

-0. 9 

0. 12 

0. 00 

7 

O 

0.3 

0 . 09 

0 . 00 

0 . 3 

0 . 09 

(!) . 00 

7 

o 

-0.0 

0. 00 

0. 00 

0. 4 

0. 10 

0 . 00 

7 

4 

0.3 

0o 10 

0 . 00 

0. 2 

0. 10 

0 . 00 

7 

5 

-0. 0 

0. 00 

0. 00 

0. 4 

0. 12 

0. 00 

7 

6 

-0 . 0 

0. 00 

0 . 00 

0. 4 

0 . 1 2 

0 . 00 

7 

7 

1.8 

0. 16 

0.00 

-0 . 0 

0 . 00 

0 . 00 

cl 

0 

0.7 

0. 10 

0. 00 




s 

1 

0 . 0 

0 . 00 

0. 00 

-0 . 6 

0 . 1 1 

0 . 00 



0.3 

0 . 09 

0 . 00 

-0 . 0 

(!) - 00 

(!) . 00 

i“i 


-0 . o 

0 . 00 

0. 00 

0 . 0 

0. 00 

0 . 00 

o 

c 

4 

-0.9 

0. 10 

0 . 00 

-0 . 0 

0 . 00 

0 . 00 

o 

5 

0.3 

0 . 03 

0 . 00 

0.7 

0 . 09 

0 . 00 

c 

h 

0 . 0 

0 . 00 

0 . 00 

-•0 . 9 

0 . OS 

0 . 00 

c 

c* 

7 

-0 . 0 

0 . 00 

0. 00 

-0 . 0 

(!) . 00 

(!) . 00 

c» 

1*1 

•3 

4M • A. 

0. 15 

0. 00 

-1 . 7 

0. 19 

(!) . 00 

9 

0 

-0 . 3 

0.11 

0 . 00 




9 

1 

0 . 0 

0 . 00 

0. 00 

0 . 5 

0 . 03 

(!) . 0< ' 

9 

2 

0.5 

0 . 03 

0 . 00 

-0 . 0 

0, 00 

0 . O' 

9 


0 . 0 

0. 00 

0 . 00 

0 . 9 

0 . 09 

0. 00 

9 

4 

0 . (I) 

0 . 00 

0 . 00 

0 . 6 

0. 07 

(!) . (DO 

9 

3 

-0. 1 

0 . 06 

0. 00 

-0 . 0 

0 . (!)(.) 

0. 00 

9 

6 

0 . 0 

0 . 00 

0 . 00 

-o . 0 

(!) . 00 

0 . 0(. 

9 

7 

-0 . 0 

0 . 00 

0. 00 

-0 . 0 

0 . 00 

0. oo 

9 

o 

-0 . 4 

0.11 

0 . 00 

0. 4 

0. 10 

0 . 00 

9 


-t . 1 

0. 16 

0. 00 

2. 1 

0. 16 

0 . (DO 


t'ri ; ♦ ~ P ' 

l.'iftii-'J ■»*■> 

OF POOR 


VttJX £$ 

quautv 
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Table 6 Con t ' d 


n 

in 

9„" 

Standard 

Daun-Ousk 


s t andard 

Oawn-Ousk 




error 

Oi t Terence 


error’ 

Oi fference 

10 

o 

0 . '9 

0.07 

0 . 00 




10 

1 

0 . 0 

(!) . 00 

0 , 00 

-0 . 3 

0 . 0 7 

0 . 00 

to 


-0. 5 

0 . 03 

0 . 00 

-0 . 0 

0 . 00 

0 . 00 

to 


-0 . .3 

0. 06 

0 . 00 

-0 . 3 

0 v 0 7 

0 . 00 

to 

4 

0. 7 

('.) . 06 

0 . 00 

-0 . 

0 . 06 

(.i . o< > 

10 

5 

0 . 0 

0 . 00 

0 . 00 

-0 . 

Ci . 03 

0 . 00 

10 

6 

0 . 0 

0 . 00 

o . oo 

-u. ll 

Co 00 

0 . oo 

to 

7 

-0. 7 

0.04 

0 . 00 

-0 . 0 

0 . 00 

0 . oo 

to 

‘3 

1.3 

0.07 

0 . 00 

-0 . 0 

0,0C) 

0 . 00 

10 

•? 

0. 1 

0.03 

0 . 00 

-0 . 3 

0 . 0* 

0 . 00 

10 

10 

-0 . 5 

0. 12 

0.00 

-0 . 0 

0.00 

0 . 00 


V 


* 4 

a 
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Table 7: Statistics of Data Sets relative to 0SFC( 12/83) 

Units are nT. 




R 

X 

Y 

l 

Dawn data 


RMS 

6.8 

8.9 

8.3 

7.3 


Mean 

.5 

.2 

-.8 

-4.6 


Si gma 

6,8 

8.9 

8.3 

5.7 

Dusk data 


RMS 

8.4 

12.8 

14.2 

o.R 


Mean 

-1.6 

-3.5 

4.3 

-.7 


Sigma 

8.3 

12.3 

13.5 

9.8 

Merged dawn and dusk 

(corrected dusk data) 

RMS 

7.0 

8.2 

8.1 

6.9 


Mean 

.5 

-1.1 

1.3 

-3.7 


Sigma 

7.0 

8.1 

8.0 

5.9 

Observatory Data 
(1978.5 - 1982.5 

with biases) 

RMS 

12.9 

16.9 

9.4 

14.4 


Mean 

-.4 

-.1 

.02 

-.2 


Sigma 

12.9 

16.9 

9.4 

14.4 

(1978.5 - 1982.5 

without biases) 

RMS 

396.7 

212.8 

211.5 

379.9 


Mean 

4.1 

3.1 

13.7 

-16.2 


Sigma 

396.7 

212.8 

211.0 

379.6 


ORIGINAL PACS P 
OF POOR QUALSVY 


Tabu 8: 

Compari son 

of g^ 0 for 

several Magsat models 
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Model 

Egoch 

9 1 °<nT) 

Difference from Reference 

GSFCC12/83) at 
same egoch (nT)_ 

MGST (6/80) 
GSFCC9/8Q) 
IGRF 1980 

1979.85 

1980.0 

1980.0 

-29989.6 

-29987.9 

-29988.0 

-.9 Lange l et a l . , 1 980 
-3.7 Langel et al., 1982 
-3.6 Peddi e, 1932 



tj? 




| CVs ! 

% 

MLftl’V 



Table 9: Statistics 

of Data 

Sets relative to 

IGRF80 

with external 

field 

parameters 

from GSFC( 12/83) . Units 

are nT. 





B 

X 

Y 

Z 

Dawn data 

RMS 

15.7 

11.7 

11.0 

13.3 


Mean 

-4.6 

- • 6 

-.7 

-4.2 


Sigma 

15. n 

11.7 

11.0 

12.6 

Dusk data 

RMS 

16. n 

15.0 

16.2 

14.9 


Mean 

-7.1 

-4.5 

4.2 

-.7 


Sigma 

14.3 

14.3 

15.7 

14.9 

Merged dawn and dusk 






(corrected dusk data) 

RMS 

14.8 

11.6 

11.0 

13.2 

Mean 

-4.1 

-2.0 

1.3 

-3.6 


Sigma 

14.2 

11.4 

10. 9 

12.7 

Observatory Data 
(1978.5 - 1982.5 






with biases) 

RMS 

26.3 

24.8 

19.9 

30.9 


Mean 

-3.8 

-5.6 

3.9 

.6 


Sigma 

26.0 

24.2 

19.5 

30.9 

(1978.5 - 1982.5 






without biases) 

RMS, 

398.4 

211.8 

214.0 

381.7 


Mean 

.6 

-2.3 

17.7 

-15.4 


Sigma 

398.4 

211.8 

213.3 

381.4 
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Table 11: 

Statistics relating 

to anomaly bias determination 


Component 

Rms of 
biases from 
GSFC(9/80), nT 

Rms of 
biases from 

GSFC( 12/83) , nT bias 

Rms of 
di fferences 

X 

292 

210 

32 

Y 

344 

199 

38 

Z 

507 

386 

50 


Figure 1 


Figure 2 


Figure 3 


Figure 4 


Figure 5 


Figure 6 


pm™ ® 

0F QUALITY 


Figure Captions 


: Difference between computed fields from the dawn model and the dusk 

model at 500 km altitude. Units are nT. (a) X-component. 

(b) Z-component. 

: Difference between computed fields from the dawn model and the 

corrected dusk model at 500 km altitude. Units are nT. 

(a) X-component. (b) Z-component. 

: Secular variation at 1980 from the GSFC( 12/83 ) model. 

(a) 5, (b) X, (c) Y, (d) 2, (e) A . Units are nT /yr. 

: Distribution of Magnetic Observatories used in the GSFC ( 12/83 ) 

solution. 

: X , Y and Z component values from observatories for 1978-1982. 

Data points are annual means; solid line is computed from 
GSFC ( 1 2/83 ) ; dashed line is computed from IGRF 19LG. The model 
components (solid and dashed lines) for X , Y and Z are 
identified by (x) , (7) , and (z) , respecti vel y . Both models use 
the observatory biases computed with GSFC (12/33). Units are nT. 

(a) Boulder, (b) Guam, (c) Huancayo , (d) Muntinlupa, (e) Tromso. 

Spherical harmonic coefficients from selected models versus time 
from 1960-1980. Units are nT. The selected models are listed in 
Table 10. (a) g* 0 , (b) g^ 0 , (c) li; 1 , (d) li, 1 . 


X-COMPONENT, nT 



DAWN-DUSK MODELS; EPOCH: 1980 
AT 500km ALTITUDE 








r igure 2a 



























Figure 5b 




HUANCAYO 

LAT -12.05 LON -75.34 ALT 3.31 Km 



ME 



Figure 5d 

MUNTINLUPA 

LAT 14.38 LON 121.02 ALT 0.06 Km 




TROMSO 

LAT 69.66 LON 18.95 ALT 0.11 Km 
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